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Diffusion of impurities in a granular gas undergoing fiomogeneous cooling state is studied. The 
results are obtained by solving the Boltzmann-Lorentz equation by means of the Chapman-Enskog 
method. In the first order in the density gradient of impurities, the diffusion coefficient D is 
determined as the solution of a linear integral equation which is approximately solved by making 
an expansion in Sonine polynomials. In this paper, we evaluate D up to the second order in 
the Sonine expansion and get explicit expressions for D in terms of the coefficients of restitution 
for the impurity-gas and gas-gas collisions as well as the ratios of mass and particle sizes. To 
check the reliability of the Sonine polynomial solution, analytical results are compared with those 
obtained from numerical solutions of the Boltzmann equation by means of the direct simulation 
Monte Carlo (DSMC) method. In the simulations, the diffusion coefficient is measured via the 
mean square displacement of impurities. The comparison between theory and simulation shows 
in general an excellent agreement, except for the cases in which the gas particles are much heavier 
and/or much larger than impurities. In theses cases, the second Sonine approximation to D improves 
significantly the qualitative predictions made from the first Sonine approximation. A discussion on 
the convergence of the Sonine polynomial expansion is also carried out. 



I. INTRODUCTION 

The use of hydrodynamic-like type equations to describe the macroscopic behavior of granular media under rapid 
flow conditions has been widely recognized in the past few years. The main difference from ordinary fluids is the 
dissipative character of collisions which leads to modifications of the Navier-Stokes equations. In order to gain 
some insight into the derivation of hydrodynamics from a kinetic theory point of view, an idealized model is usually 
considered: a system composed by smooth hard spheres with inelastic collisions. In the low-density regime, the 
above issue can be addressed by using the Boltzmann kinetic theory conveniently modified to account for inelastic 
collisions 0. In this case, assuming the existence of a normal solution for sufficiently long space and time scales, 
the Chapman-Enskog method can be applied to solve the Boltzmann equation and determine the explicit form 
of the transport coefficients. This objective has been widely covered in the case of a monocomponent gas |^^|^, 
where all the particles have the same mass and size. Nevertheless, a real granular system is generally composed by 
particles with different mass densities and sizes, which leads to many interesting phenomena observed in nature and 
experiments such as separation or segregation. 

Needless to say, the determination of the transport coefficients of a multicomponent granular fluid is much more 
complicated than in the case of a monocomponent system. Several attempts to get those coefficients from the 
Boltzmann equation began time ago IflHIMSIlfl i but the technical difficulties for the analysis entailed approximations 
that limited their range of validity (quasielastic limit). In particular, all these works assume energy equipartition 
and consequently, the partial temperatures Ti are made equal to the global granular temperature T. However, the 
failure of ei iergy equip artition in granular fluids 0, 0, 0, 0, has also been confirmed by computer simulations 
Ir^ ITsL Il9l I20I I21I |23| and even observed in recent experiments |230|. As a consequence, previous studies for 
granular mixtures based on a single temperature must be reexamined by theories which take into account the effect 
of temperature differences on the transport coefficients. For this reason, recently Garzo and Dufty '2^1 have carried 
out a derivation of hydrodynamics for a granular binary mixture at low-density that accounts for nonequipartition 
of granular energy. Their results provide a description of hydrodynamics in granular mixtures valid a priori over 
the broadest parameter range and not limited to nearly elastic particles. On the other hand, from a practical point 
of view, the expressions for the transport coefficients were obtained by considering the leading terms in a Sonine 
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polynomial expansion of the distribution functions. In the case of the shear viscosity coefScient, the first Sonine 
predictions compare quite well with a numerical solution of the Boltzmann equation in the uniform shear flow state 
|2(t| | obtained from the direct simulation Monte Carlo (DSMC) method [23| ■ Exceptions to this agreement are extreme 
mass or size ratios and strong dissipation. These discrepancies could be mitigated in part if one considers higher-order 
terms in the Sonine polynomial expansion. 

The evaluation of the transport coefficients for a granular binary mixture beyond the first Sonine approximation is 
quite a hard task, due mainly to the coupling among the different integral equations associated with the transport 
coefficients (25.| . For this reason, to make some progress into this problem, one needs to study some specific simple 
situations. Here, we consider the special case in which one of the components (say, for instance, species 1) is present 
in tracer concentration. The tracer problem is more amenable to be treated analytically since the tracer particles 
(impurities) are enslaved to the granular gas and there are fewer parameters. In this situation, one can assume that 
the granular gas (excess component) is not affected by the presence of impurities so that its velocity distribution 
function /2 verifies the (closed) nonlinear Boltzmann equation. Moreover, the mole fraction of impurities is so small 
that one can also neglect collisions among tracer particles in their corresponding kinetic equation. As a consequence, 
the velocity distribution function /i of impurities verifies the linear Boltzmann-Lorentz equation. 

The goal of this paper is to determine the diffusion coefficient of impurities immersed in a low-density granular 
gas undergoing the so-called homogeneous cooling state (HCS). It corresponds to a homogeneous state where the 
temperature uniformly decreases in time due to dissipation in collisions. Diffusion of impurities in a homogeneous gas 
is perhaps the simplest nonequilibrium problem one can think of. The diffusion coefhcient D is obtained by solving the 
Boltzmann-Lorentz equation by means of the Chapman-Enskog method through the first order in the concentration 
gradient. As in the elastic case, the coefficient D is expressed in terms of the solution of a linear integral equation 
that can be solved by making an expansion in Sonine polynomials. An important simplification with respect to the 
case of arbitrary composition 25] is that the diffusion coefficient satisfies a closed equation and so is not coupled 
to the remaining transport coefficients. This simplifies enormously the procedure of getting the different Sonine 
approximations. Here, we have retained up to the second Sonine approximation and have determined D in terms of 
the coefficients of restitution for the impurity-gas and gas-gas collisions, and the parameters of the system (masses 
and sizes). As for the shear viscosity coefficient 5^, kinetic theory predictions are also compared with numerical 
solutions of the Boltzmann-Lorentz equation by using the DSMC method. In the simulations, the diffusion coefficient 
is computed from the mean-square displacement of impurities once a transformation to dimensionless variables allows 
one to get a stationary diffusion equation. 

Some relat ed p apers studying diffusion in a homogeneous granular gas have been published earlier. Thus, Brilliantov 
and Poschel [2^ obtained the self-diffusion coefficient (when impurities and gas particles are mechanically equivalent 
particles) from a cumulant expansion of the velocity autocorrelation function, while Brey et al. 29] evaluated this 
coefficient in the first Sonine approximation by means of the Chapman-Enskog expansion. In the latter work, the 
authors also performed Monte Carlo simulations and found a good agreement between theory and simulation. The 
dynamics of a heavy impurity particle in a gas of much lighter particles has been also studied from the Fokker-Planck 
equation '30]. The predictions of this (asymptotic) theory have also been confirmed by Monte Carlo and molecular 
dynamics simulations 31]. More recently, Dufty et al. p33 ] applied the familiar methods of linear response theory 
to get Green-Kubo and Einstein representations of diffusion coefficient in terms of the velocity and mean-square 
displacement correlation functions. These correlation functions were also evaluated approximately considering the 
leading cumulant approximation. This result coincides with the one obtained from the Chapman-Enskog method in 
the first Sonine approximation j^. The expression for the diffusion coefficient obtained in Ref. [HI has been also 
compared with results from molecular-dynamics simulations |34( over a wide range of density and inelasticity, for 
the particular case of self-diffusion. The comparison shows that the approximate theory is in good agreement with 
simulation data up to moderate densities and degrees of inelasticity, although large discrepancies exist at high density 
and large dissipation. However, these deviations from the Enskog model are primarily due to the effect of an unstable 
long wavelength shear mode in the system. Our work complements and extends the results obtained in previous 
studies in two aspects; First, our Chapman-Enskog solution is more accurate since incorporates higher order Sonine 
polynomials, and second, our simulations cover a range of values of the ratios of mass and particles sizes where the 
effect of energy nonequipartition on diffusion is in general quite important |25| . The diffusion phenomenon in granular 
shear flows has also been widely studied by computer simulations , kinetic theory Isij and real experiments (37| . 
In this situation, due to the presence of shear flow, the resulting diffusion process is anisotropic and, thus, it must be 
described by a diffusion tensor. 

The plan of the paper is as follows. In Sec^Jwe describe the problem we are interested in and analyze the state of the 
granular gas and impurities in the absence of diffusion. In particular, we show the breakdown of energy equipartition 
and illustrate the dependence of the temperature ratio on the parameters of the problem. Section IlIII deals with the 
Chapman-Enskog method to solve the Boltzmann-Lorentz equation in the first order of the concentration gradient. 
Some technical details of the calculations are given in Appendices IbI and lO In Sec. El we present the Monte Carlo 
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simulation of the Boltzmann-Lorentz equation and compare the simulation data with the theoretical results obtained 
in the first and second Sonine approximations. Finally, in Sec. [V]we close the paper with a brief discussion on the 
results obtained in this paper. 

II. GRANULAR BINARY MIXTURE IN THE HOMOGENEOUS COOLING STATE 

Consider a binary mixture of smooth hard spheres of masses mi and TO2, diameters ai and CT2, and interparticle 
coefficients of restitution an, 022, and — ck2i- Here, a^j is the coefficient of restitution for collisions between 
particles of species i and j. In the low-density regime, the distribution functions /i(r,v;<) (i = 1,2) for the two 
species are determined from the set of nonlinear Boltzmann equations 1] 

{dt + VI . V) /. = ^ J,, [vi|/.(t), f,{t)] , (1) 
j 

where the Boltzmann collision operator Jij [vi|/i, fj] is 

JijWilftJj] = CTtj J j da Q{a ■ gi2){a ■ ^12) 

X Kr2/^(r,v'i,t)/,(r,v^,i)-/,(r,vi,t)/,(r,V2,i)] . (2) 

Here, = (cr; + Uj) /2, is a unit vector along their line of centers, is the Heaviside step function, and gi2 = 
Vi — V2. The primes on the velocities denote the initial values {v'i,V2} that lead to {vi,V2} following a binary 
collision: 

v'l = vi - (1 + a^r.i) (ct . gi2)CT, V2 = V2 + (1 + %-^) (5^ • gi2)5^, (3) 
where /^y = m^/ (m^ + rrij). 

The relevant hydrodynamic fields are the number densities Ui , the flow velocity u, and the "granular" temperature 
T . They are defined in terms of moments of the distributions /j as 

n^ = j d-vf,{-v) , pu = Y^ J dvm.,vf.,{-v) , (4) 

nT = p = J2j dv^V^V.(v), (5) 

where n = ni + 712 is the total number density, p = mini + m2Ti2 is the total mass density, p is the hydrostatic 
pressure, and V = v — u is the peculiar velocity. Furthermore, it is convenient to introduce the kinetic temperatures 
Tj for each species, which measure their mean kinetic energies. They are defined as 

^n,r, = J dv^^V., (6) 

so that the granular temperature T is 

T = ^x,r„ (7) 

i 

where Xi = rii/n is the mole fraction of species i. 

The collision operators conserve the particle number of each species and the total momentum but the total energy 
is not conserved: 

J d^rJ^j [v|/„ /,] =0, Y.J o'vm.v J,, [v|/„ /,] = 0, (8) 
3nTC - - E / rfvm,T/V„ [v|/„ f,], (9) 
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where C is identified as the cooling rate due to inelastic colhsions among all species. At a kinetic level, it is convenient 
to characterize energy transfer in terms of the cooling rates associated with the partial temperatures T^. They are 
defined as 

5n^T,Q = -Y. [ c^vm,FV,,[v|/„ /,]. (10) 
The total cooling rate C, can be expressed in terms of the partial cooling rates C,i as 

C = r-i^a:,r,c.. (11) 

i 

Because of the complexity embodied in the general description of a binary mixture, here we consider the special 
case in which the mole fraction of one of the components (say, for instance, 1) is negligible [xi <C 1). We are interested 
in studying the diffusion of impurities moving in a background granular gas undergoing homogeneous cooling state 
(HCS). In the tracer limit, one expects that the state of the granular gas (solvent) is not disturbed by the presence 
of impurities (solute) and so in all of the following it is assumed that the gas is in its HCS. In addition, collisions 
among impurities themselves can be neglected versus the impurity-gas collisions. Under these conditions, the velocity 
distribution function /2 of the gas verifies a (closed) Boltzmann equation and the velocity distribution function /i of 
impurities obeys a (linear) Boltzmann-Lorentz equation. Let us start by describing the state of the system in the 
absence of diffusion. 



A. Granular gas 



As said above, the gas is in the HCS. This state corresponds to a homogeneous solution of the nonlinear Boltzmann 
equation, /2(w,t), in which all the time dependence occurs through the temperature of the gas T2{t) ~ T{t). In this 
case, the time derivative of /2 can be represented more usefully as 



dth ^ ~(:2TdTh = \q2^ ■ i^h) , 

where use has been made of the balance equation for the temperature 

T-^dtT = -C2, 

with 

1 



C2 = - 



3^2?" 



d-Vm2V J22[./2,/2] 



The Boltzmann equation can be written as 



1 d 



(12) 
(13) 
(14) 

(15) 



The second equality in H12|) follows from dimensional analysis which requires that the temperature dependence of /2 
occurs only through the temperature T{t). Consequently, f2{v,t) has the form 



f2{v,t) = n27r-^^^v-''/\t)<i>2{v/vo{t)), 



(16) 



where vo{t) — yJ^XfJirvi is the thermal velocity of the gas. So far, the exact form of $2 has not been found, although 
a good approximation for thermal velocities can be obtained from an expansion in Soninc polynomials. In the leading 
order, $2 is given by 



$2(W*) 



*4 



. *2 15 
5w + -r 



(17) 



The coefficient C2 can be obtained by substituting first Eqs. H16|) and H17|) into the Boltzmann equation (|15|1 . multiplying 
that equation by and then integrating over the velocity. When only linear terms in C2 are retained, the estimated 
value of C2 is 381 



02(022) = 



32(l-a22)(l-2ai2) 
81- 17a22 + 30ai2(l -"22)' 



(18) 
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FIG. 1: Temperature ratio 7 — Ti/T versus the coefficient of restitution a = ai2 — 0:22 for /i = m\/m2 = 2 and uj = (T\/a2 ~ 1 
(open circles) and /i = 1/4 and uj — 1/2 (filled circles). The lines are the theoretical results and the symbols refer to the 
numerical results obtained from the DSMC method. 



Additionally, the cooling rate C2 can also be determined from the Sonine approximation H17|) . The result is [36 

C2 = -\/27rn2cr2«o(l - "22) ( 1 + ^C2 
Estimate (|18|l presents quite a good agreement with Monte Carlo simulations of the Boltzmann equation [39l |4 



(19) 



B. Impurities 



In the absence of diffusion, impurities are also in HCS and so its velocity distribution function fi(v,t) satisfies the 
Boltzmann-Lorentz equation 



The time evolution of Ti (t) is 



where 



^Cl|;-(v/l) = Jl2[/l,/2] 



T^^dtTi = -Ci 



Ci 



dvmiv Ji2[/i,/2] 



From Eqs. (|13|l and H21|l one easily gets the time evolution of the temperature ratio 7 = Ti/T: 

I'^da = C2 - Ci- 



(20) 



(21) 



(22) 



(23) 



The fact that /i depends on time only through T{t) necessarily implies that the temperature ratio 7 must be inde- 
pendent on time, and so Eq. (|23|) gives the HCS condition C,i{t) — C,2{t) = C(0- However, although both components 
have a common cooling rate, their partial temperatures are different. In other words, the impurity equilibrates to a 
common HCS with different temperatures for the impurity and gas particles. This implies a breakdown of the energy 
equipartition. The violation of energy equipartition in multicomponent granular systems has been even observed in 
real experiments of vibrated mixtures in two 23] and three 24] dimensions. 

As for the gas, dimensional analysis requires that the solution to Eq. H20|l is of the form 



(24) 



The determination of fi to leading order in the Sonine expansion has been analyzed elsewhere for arbitrary composition 
|llj ] and only the main results are quoted here. In the first Sonine approximation, the distribution $1 is given by 



$i(w*) ^0^/2 



^ ( - Wv*^ + - 
4 V 4 



(25) 
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where 9 = miT/m2Ti is the mean-square velocity of the gas particles relative to that of the tracer particles. The 
coefficient Ci in Eq. H25|) can be determined by substitution of Eqs. (|16|) and (I24|) into Eq. H2U|I and retaining all terms 
linear in ci and C2 for the leading polynomial approximations 1)17(1 and ((25|l . Once the coefficient ci is known, one 
can estimate the temperature ratio 7 = Ti/T from the constraint C,i = ^2- The solution to this equation gives 7 as a 
function of the mass ratio /i = mi/m2, the size ratio w = cti/ct2 and the coefficients of restitution 022 and ai2. The 
explicit expressions of the partial cooling rate C,i and the coefficient ci are displayed in Appendix IXI 

Except for some limiting cases (elastic case and mechanically equivalent particles), our results yield 7 7^ 1, and 
so the total energy is not equally distributed between both species. The lack of energy equipartition has dramatic 
consequences in the large impurity/gas mass ratio since there is a peculiar "phase transition" for which the diffusion 
coefficient is normal in one phase and grows without bound in the other [4ll | . To illustrate the violation of equipartition 
theorem, in Fig. ^ we plot the temperature ratio versus the coefficient of restitution a for two different cases: /i = 1/4, 
LD ~ 1/2, and /i = 2, a; = 1. For the sake of simplicity, we have taken a common coefficient of restitution, i.e., 
a = ai2 — a22- We also include the simulation data obtained via a numerical solution of the Boltzmann equation by 
means of the direct simulation Monte Carlo (DSMC) method ■ It is apparent the excellent agreement found between 
the theory and simulation, showing the accuracy of the approximations 1(17(1 and 1(25(1 to estimate the temperature ratio. 
We also observe that the temperature of the impurity is larger than that of the gas when the impurity is heavier than 
the grains of the gas. In addition, the deviations from the energy equipartition increase as the mechanical differences 
between impurities and particles of the gas increase. As we will show later, in general the effect of the temperature 
differences on the diffusion coefficient is quite important. 



III. DIFFUSION COEFFICIENT 



We want to determine the diffusion coefficient of tracer particles immersed in a granular gas in HCS. The diffusion 
process is induced by a weak concentration gradient Vxi, which is the only gradient present in the system. Under 
these conditions, the kinetic equation for /i reads 

(5t+vi-V)/i = Ji2[vi|/i,/2] . (26) 

Impurities may freely exchange momentum and energy with the particles of the granular gas, and, therefore, these 
are not invariants of the Boltzmann-Lorentz collision operator J12 [fi , /2] . Only the number density of impurities is 
conserved: 

dtni + = 0, (27) 

mi 

where the flux of tracer particles ji is deflned as 

ji = mi / dvv/i(v). (28) 



The conservation equation 1(27(1 becomes a closed hydrodynamic equation for ni once ji is expressed as a functional 
of the fields ni and T. Our aim is to get the mass flux in the flrst order in Wxi by applying the Chapman-Enskog 
method The Chapman-Enskog method assumes the existence of a normal solution in which all the space and 
time dependence of fi occurs through the hydrodynamic flelds. In the tracer diffusion problem, this means that 

/i(r,v;t) =/i[v|.Ti(r,i),T(t)], (29) 

where it has been assumed that /2 also adopts the normal form ((16() . The Chapman-Enskog procedure generates the 
normal solution explicitly by means of an expansion in gradients of the fields 

h = A'^+efi'^+e'fi'^ + ---, (30) 

where e is a formal parameter measuring the nommiformity of the system. Given that in our problem the only 
inhomogeneity is in the mole fraction xi, each factor e corresponds to an implicit Vxi factor. The time derivative is 
also expanded as dt — d^^^ + edj:^^ + • • • , where 

di°'^xi = 0, V = -TC, (31) 



■J 



(fe-i) 



TOirt2 



dl'^^T = 0, A: > 1, 



(32) 
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with 



j[-> =m,J dvv/r- (33) 

Upon deriving Eqs. H31|l and (|32|l . use has been made of the balance equations and H27(l and the constraint 

(^f) = C2 = C- Here, cf '> is given by Eq. ^ with A ^ f^°\ 

The zeroth-order approximation /{°'' is the solution of (|20|l . whose approximate form is given by Eqs. (|24|l and l|25|) 

but taking into account now the local dependence on the mole fraction xi. Since is isotropic, it follows that the 

flux of impurities vanishes at this order, i.e., j^^'' = 0, and so d^^^xi = 0. To first order in e, one has the kinetic 
equation 



(0) 



A;<».)v,.v.. ,34, 

The second equality follows from the balance equations and the space dependence of /}°^ through xi. The 
solution to Eq. H34|) is proportional to Vxi, namely, it has the form 

f[^^^A-Wxi. (35) 

The coefficient ^ is a function of the velocity and the hydrodynamic fields xi and T. According to (|31|l . the time 

derivative cJj"'' acting on A. can be evaluated by the replacement —QTOt where Q is given by (|19|l . Thus, 

substitution of into yields 

-CTBtA ~ J12 [A, h] = - [J^/i^) V- (36) 
To first order, the mass flux has the structure |42| 

= -miD^xi, (37) 

where D is the diffusion coefficient. Use of Eq. H35() into the definition (|33() (for fc = 1) allows one to identify the 
coefficient D. It is given by 

D = -^J dvvA. (38) 

For practical purposes, the linear integral equation (|3()|l can be solved by using a Sonine polynomial expansion. Our 
goal here is to determine the diffusion coefficient up to the second Sonine approximation. In this case, the quantity 
A is approximated by 

A^/i,M[aiv + a2Si(v)], (39) 
where /i,m is a Maxwellian distribution at the temperature Ti of the impurities, i.e. 

/»/(v)^n,feV"»p(-=^V (40) 

and Si(v) is the polynomial 



Si(v)== ( imi«2_ (41) 



The coefficients ai and 02 are defined as 



a^ = ^ fd...A^-^, (42) 
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FIG. 2; Plot of the function A = _D[2]/_D[1] versus the coefficient of restitution a for three binary mixtures: (a) \i — 1/8 and 
Lo = 1/2; (b) = 125/64 and lu = 5/4; and (c) = 8 and u = 2. 
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(43) 



These coefficients are determined by substitution of Eq. H39I) into tlie integral equation (|36|l . The details are carried 
out in Appendix IbI The second Sonine approximation -D[2] to the diffusion coefficient can be written as 



D\2\ = D\\\- 



(^^:-K*)[^^d-|r-(ci/27K] 



{K-\^*){^l-\C)-K\K-{^*h)\ 



^ i?[l]A, 



(44) 



where the function A can be easily identified from the second equality and -D[l] refers to the first Sonine approximation 
to the diffusion coefficient. Its expression is 



D{\\ 



1 



(45) 



Here, vq = nia^vo, 



Vq 6 



— Co 
32 



(46) 



and the quan tities , v'^^ v*, and given in Appendix [CI Expression (|45|l coincides with the one previously 

obtained [3^ |^ from a Green-Kubo formula in the leading order in a cumulant expansion. 

In general, D[l] and D[2] present a complex dependence on the coefficients of restitution and the ratios of mass 
and sizes. Before analyzing this dependence, it is instructive to consider some special cases. In the elastic limit 
(ck22 — Oi\2 = 1), one has ci = C2 = 0, 7 = 1, and the first i?o[l] and second -Do[2] Sonine approximations to the 
diffusion coefficient become 



3 ^(mi +m2)T 
-t-'oilj = -r- 



(47) 



^ ^ 30/i2 + 16/i + 12 



(48) 



Equations (|47|l and H48|) coincide with the ones derived by Mason for a gas- mixture of elastic hard spheres. 
Moreover, in the case of mechanically equivalent particles (mi = ra-i^ (j\ = <J2, 0,22 = 0^12 = a), Eq. (|45() reduces to 
the one derived in the self-diffusion problem 29]. 

The function A = Z?[2]/Z3[l] is plotted in Fig. |21 versus the coefficient of restitution a for three different values of 
the mass ratio /i and the size ratio lo. We consider a binary mixture where the mass density of the impurity is equal 
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to that of a bath particle and so, /z = w^. As in the elastic case, we see that A clearly differs from 1 when impurities 
are much lighter than the particles of the gas. This means that the widely used first Sonine approximation is not 
sufficiently accurate for this range of values of mass and diameter ratios. On the other hand, the convergence of the 
Sonine polynomial expansion improves when increasing the mass and size ratios, and the first Sonine approximation 
is expected to be quite close to the exact value of the diffusion coefficient. These conclusions arc qualitatively similar 
to those found in the case of elastic hard-sphere mixtures 44!| . 

To close this section, let us write the diffusion equation. In the hydrodynamic regime (where the normal solution to 
the Boltzmann-Lorentz equation holds), the diffusion coefficient D{t) depends on time only through its dependence 
on the temperature T{t). According to Eqs. and D{t) oc ^jT{t). This time dependence can be eliminated 
by introducing appropriate dimensionless variables. A convenient set of dimcnsionless time and space variables is 
given by 

r= fdt'Mt'), (49) 

Jo Mt) 

The dimensionless time scale r is the integral of the average collision frequency and thus is a measure of the average 
number of collisions per gas particle in the time interval between and t. The unit length Vf){t) / vo{t) introduced in 
the second equality of H49|) is proportional to the time- independent mean free path of gas particles. In terms of the 
above variables, the diffusion equation l|?7)l becomes 

^-^=D*Vlx„ (50) 

OT 



where is the Laplace operator in I space and 



^ (51) 



While D{t) (X y/T{t), the reduced diffusion coefficient D* is time- independent. This is closely related with the validity 
of a hydrodynamic description of the system for times large compared with the mean free time. In this context, Eq. 
(|50|l is a diffusion equation with a constant diffusion coefficient D* . It follows that the mean square deviation of the 
position £ of impurities after a time interval r is j42j 

{\e{T)-m\^)^6D*r. (52) 

Restoring the dimensions to the average position and time, Eq. H52|l can be rewritten as 

Equation H53|) is the Einstein form, relating the diffusion coefficient to the mean-square displacement. This relationship 
will be used later to measure the diffusion coefficient by means of the DSMC method. 



IV. COMPARISON BETWEEN THEORY AND MONTE CARLO SIMULATIONS 



The adaptation of DSMC method to analyze binary granular mixtures has been described in detail in previous 
works (see, for instance, Ref. ^3), so that here we only shall mention the aspects related to the specific problem of 
diffusion of impurities in a granular gas under HCS. The main difference between the procedure used here and the one 
described in Ref. ,16.] is its restriction to the tracer limit {xi — > 0). Because of this, during our simulations collisions 
1-1 are not considered, and when a collision 1-2 takes place, the postcollisional velocity obtained from the scattering 
rule is only assigned to the particle 1. According to this scheme, the numbers of particles Ni have simply a statistical 
meaning, and hence they can be chosen arbitrarily. 

In the course of the simulation one can distinguish two stages. In the first one, the system (impurities and gas 
particles) evolves from the equilibrium initial state to the HCS. In the second stage, the system is assumed to be 
in the HCS, and then the kinetic temperatures Ti{t) and the diffusion coefficient D{t) are measured. The latter is 
obtained from the mean square deviation of the position of impurities [Eq. H53|) ]. i.e., 



(54) 
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FIG. 3: Plot of D/Do[l] as a function of the dimensionless time r for ^ — 2, uj — 1 and a = 0.7, 0.8 and 0.9. Here, Do[l] refers 
to the first Sonine approximation to the elastic value of the diffusion coefficient. 



Here, \ri{t) — ri(0)| is the distance traveled by the impurity i from t = until time t, t ~ being the beginning of the 
second stage. In addition, (• ■ ■ ) denotes the average over the iVi impurities and At is the time step. In our simulations 
we have typically taken 10^ particles of the granular gas, 5 x lO"' impurities, and a time step At = 3 x 10~^A/wo, 
where A = {^/2Trn2(J2)''^ is the mean free path for collisions among granular gas particles. The results are averaged 
over a number of Af — b replicas. 

If the hydrodynamic description (or normal solution in the context of the Chapman-Enskog method) applies one 
expects that, after a transient regime in which each gas particle has collided about five times, the reduced diffusion 
coefficient 



D{t) ^ f<Ji2^ ' 




(55) 



achieves a time-independent plateau. This is illustrated in Fig. 13 for a system with /i = 2, w = 1 and three values of 
a : a = 0.7, 0.8 and 0.9. The dimensionless time scale r is defined in Eq. (|49|l and measures the average number of 
collisions per gas particle, which is of the same order as the corresponding value for impurities. We observe that after 
several collisions, D/Dq[1\ reaches a stationary value whose time average is the simulation result for the (reduced) 
diffusion coefficient. 

The steady state values of D/ Da[l] obtained from simulation data (which are calculated by averaging over a time 
period) can be compared with the theoretical predictions D[2]/£)o[l]j Eq. and Z)[l]/£'o[l], Eq. for different 
values of the parameters of the system. As in the previous figures, we assume that 0:22 = 0112 = a so that we reduce 
the parameter set of the problem to three quantities: {a,/z,iu;}. 

Let us consider first the self-diffusion case (/i = w = 1). As said above, this special situation was studied by Brey 
et al. p9l | analytically (up to the first Sonine approximation) and by computer simulations. Figure 0] shows the 
reduced diffusion coefficient D{a) / Di^[l\ as a function of the coefficient of restitution a as given by the first Sonine 
approximation (dashed line), the second Sonine approximation (solid line) and Monte Carlo simulations (symbols). 
Although the agreement between the first Sonine approximation and simulation is quite good beyond the quasielastic 
limit (say for instance, a > 0.8), discrepancies between both results increase as the value of a decreases. This 
tendency was also observed in the comparison made in Ref. poj . We observe that the agreement between theory and 
simulation improves over the whole range of values of a considered when one takes the second Sonine approximation. 
For instance, for a = 0.5 the first and second Sonine approximations to D differ by 4.7% and 1.2%, respectively, from 
the value measured in the simulation. We have also included (dotted line) the result obtained for the second Sonine 
approximation to D{a)/ Dq[1] if the distribution functions of the granular gas /2 and impurities /[^"^ at zeroth-order 
were both approximated by Gaussians, i.e., when one formally puts ci = C2 = in Eq. 144|l . We observe that the 
inffuence of these Sonine contributions to the self-diffusion coefficient is negligible, except for quite large values of 
dissipation. 

Consider now the situation in which impurities and particles of the gas can differ in size and mass. First, in Fig. 
we analyze the effect of the mass ratio on the accuracy of the Sonine approximations to the reduced diffusion 
coefficient D{a)/ Do[l]. Here, we take cti — a2 and a = 0.7. We find that, for fixed cu and a, the second Sonine 
approximation D[2\ differs from the first Sonine approximation D[l] as fi is varied. For the case considered in Fig. 
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0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
a 

FIG. 4: Plot of the reduced self-difFusion coefficient iI>/£)o[l] as a function of the coefficient of restitution a as given by the 
first Sonine approximation (dashed fine), the second Sonine approximation (soHd fine), and Monte Carlo simulations (symbols). 
Here, -Do[l] refers to the first Sonine approximation to the elastic value of the self-difFusion coefficient. The dotted line 
corresponds to the results obtained for the second Sonine approximation when one takes c\ = ci = 0. 




FIG. 5: Plot of the reduced diffusion coefficient D jDn\\\ as a function of the mass ratio /i for a = 0.7 and uj = \. The dashed 
line refers to the first Sonine approximation, the solid line corresponds to the second Sonine approximation while the symbols 
are the results obtained from Monte Carlo simulations. Here, iJofl] refers to the first Sonine approximation to the elastic value 
of the diffusion coefficient. 




0.30 0.45 0.60 0.75 0.90 



FIG. 6: Plot of the reduced diffusion coefficient -D/Do[l] as a function of the size ratio w for a = 0.7 and = 1/4. The dashed 
line refers to the first Sonine approximation, the solid line corresponds to the second Sonine approximation while the symbols 
are the results obtained from Monte Carlo simulations. Here, Do [1] refers to the first Sonine approximation to the elastic value 
of the diffusion coefficient. 
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a 



FIG. 7: Plot of the reduced diffusion coefficient D{a) / D{1) as a function of the coefficient of restitution a for jj, — 8 and 
uj — 2. The dashed line refers to the first Sonine approximation, the solid line corresponds to the second Sonine approximation 
while the symbols are the results obtained from Monte Carlo simulations. The dotted line is the first Sonine approximation 
by assuming the equality of the partial temperatures 7 = Ti/T — 1. Here, D{1) refers to the elastic value of the diffusion 
coefficient consistently obtained in each approximation. 




0.5 0.6 0.7 0.8 0.9 1.0 

a 

FIG. 8: Plot of the reduced diffusion coefficient D{a) / D{1) as a function of the coefficient of restitution a for fi = 1/8 and 
UJ = 1/2. The dashed line refers to the first Sonine approximation, the solid line corresponds to the second Sonine approximation 
while the symbols are the results obtained from Monte Carlo simulations. Here, D{1) refers to the elastic value of the diffusion 
coefficient consistently obtained in each approximation. 

|31 we see that D[2] > D[l] if /i < 0.45 while D[l] > D[2] otherwise. The comparison with simulation data shows 
again that the theoretical predictions are improved when one considers the second Sonine approximation, showing 
an excellent agreement between D[2\ and Monte Carlo simulations in the range of values of fi explored. Similar 
conclusions are found with respect to the variations of D{a)/ Dq[1] on the size ratio co at fixed values of fi and a, as 
shown in Fig. for a = 0.7 and = 1/4. Differences between both Sonine approximations are especially important 
for small values of w. In this region, the discrepancies between D[2\ and simulation data are slightly larger than the 
ones found before in Fig. [S] for the mass ratio. This means that perhaps the most significant variations of the ratio 
£'[2]/£'[l] occur when the size ratio uj is varied at fixed /x. 

Finally, we explore the influence of dissipation on the two first Sonine approximations at fixed values of /i and to. 
According to Fig.|21 it is apparent that the deviations of A = 13 [2] //^[l] from unity are important for small values of 
the mass ratio and/or the size ratio. In fact, for elastic mixtures {a = 1), A ~ 1.069 for the case {/i — 1/8, a; = 1/2}, 
while A ~ 1.00049 for the symmetric case — 8,uj — 2}. This clearly shows that in the elastic case the second 
Sonine approximation significantly differs from the first one for small values of fi and lu (Lorentz gas limit), while 
both approximations practically coincide when fi and lo are much larger than 1 (Rayleigh gas limit). In the former 
case, where the tracer diffusion coefficient can be obtained exactly, the first and second Sonine approximations to D 
differ by 13% and 5%, respectively, from the exact value of the diffusion coefficient [ij. Here, we want to analyze 
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the influence of the inelasticity on the trends already observed in ordinary fluids. In order to assess this effect, it 
is convenient to reduce the diffusion coefficient D{a) with respect to its elastic value D{1) consistently obtained 
in each approximation. Thus, D{1) Do[i-], Eq. I|47(l . when one considers the first Sonine approximation, while 
D{1) Do[2], Eq. H48|) . when one considers the second Sonine approximation. Simulation data will be also reduced 
with respect to the elastic value measured in the simulation. Figure [3 shows D{a) / D{1) for fj. — 8 and uj — 2 while 
Fig. IHlshows D{a) / D{1) for /i = 1/8 and u! = 1/2. The first finding is that the conclusions obtained in the elastic 
case on the convergence of Sonine polynomial expansion are kept at a qualitative level for granular gases: the Sonine 
polynomial expansion exhibits a poor convergence for sufficiently small values of the mass ratio ^ and/or the size 
ratio to, while this convergence improves significantly as /z and/or uj increases. As a matter of fact. Fig. |7| shows 
that both Sonine approximations are practically indistinguishable and present an excellent agreement with simulation 
data. This is not so in the case of Fig. |S1 where only the second Sonine approximation to D exhibits a qualitative 
good agreement with simulation results. As in the elastic case one perhaps would have to consider the third 
Sonine approximation to provide an accurate estimate for the diffusion coefficient. 



V. DISCUSSION 




In this paper we have analyzed the diffusion of impurities in a dilute granular gas undergoing homogenous cooling 
state (HCS). This is the simplest example of transport in a multicomponent granular gas, since the gas is in a 
homogeneous state while diffusion appears by the presence of a weak concentration gradient, which is the only 
gradient present in the system. Here, the diffusion coefficient D has been computed by combining two complementary 
approaches: a Chapman-Enskog solution to the Boltzmann-Lorentz equation and numerical solutions of the same 
equation by means of the direct simulation Monte Carlo (DSMC) method. As in the elastic case, the diffusion 
coefficient D is given in terms of the solution of an integral equation, Eq. (|3t)|) . A practical evaluation of D is possible 
by using a Sonine polynomial expansion and approximate results are not limited to weak inelasticity. Here, we have 
determined D in the first (one Sonine polynomial) and second (two Sonine polynomials) Sonine approximation as a 
function of the temperature and the different mechanical parameters of the system, namely, the (constant) coefficients 
of restitution for the impurity-gas and gas-gas collisions, the masses and the particle sizes. Our study complements 
and extends previous works on diffusion in undriven granular gases in the cases of self-diffusion |29l l33 | and Brownian 
and provides an explicit expression of the coefficient D beyond the first Sonine approximation 

|32' 

Comparison between simulation data and theory shows that in general the second Sonine approximation D[2] 
improves significantly the predictions of the first Sonine approximation D[l], especially for values of the mass ratio 
II and/or the size ratio iv smaller than 1. For this range of values of fi and uj (see, for instance. Fig. (H)) one should 
consider higher-order polynomial terms to get a quantitative good agreement with Monte Carlo simulations. However, 
the convergence of the Sonine polynomial expansion (see, for instance, Fig. |7|) improves with increasing values of fi 
and Lu and the first Sonine correction seems to be quite close to the exact value. These trends are quite similar to 
those previously found for ordinary fluid mixtures |44| . 

Apart from extreme mass or size ratios, our results show again that the accuracy of the Sonine polynomial solution 
for granular systems is similar to that for elastic collisions. Possible discrepancies between theory and simulation can 
be removed by retaining higher-order Sonine corrections. On the other hand, as in all previous studies 26, 29, 31,^, 
the good agreement found here over quite a wide range of values of the coefficient of restitution and mass and size 
ratios is an additional evidence of the validity of the hydrodynamic description to analyze some states of granular 
fluids. 

One of the main features of our theory is that it incorporates the effect of e nergy n onequipartition on diffusion. 
This new effect had already been considered in some previous studies [2^ |2^ |33, 133l33l | . To assess this effect on the 
diffusion coefficient D, in Fig. [7|we include for comparison the result for D[l] with = 8 and u) = 2 that would be 
obtained if the differences in the partial temperatures were neglected (Ti = T) . Clearly, inclusion of this effect makes 
a significant difference over the whole range of dissipation considered (for instance, the actual value of 7 is 7 ~ 2.35 
for a = 0.8). 

The results derived here are restricted to the low-density regime, small gradient in the concentration gradient, and 
to a system in the HCS. The HCS is known to be unstable under long wavelength perturbations or fiuctuations [45l |. 
with the critical wavelength Ac = 27r-^/ry72C, where 77 is the shear viscosity of the granular gas. Thus, our expression for 
the diffusion coefficient is only meaningful if the time for instability is longer than the few collision times required to 
reach the hydrodynamic regime. At low density this is clearly the case, even at strong dissipation, as demonstrated by 
both Monte Carlo and molecular dynamics simulations [2^. With respect to the extension to dense gases, the Enskog 
equation provides a useful generalization of the Boltzmann equation to higher densities for a gas of hard spheres. 
Applied to the diffusion of impurities in a granular fluid in the HCS, the only difference between the Boltzmann 
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FIG. 9: Plot of the reduced diffusion coefficient D{a) / D{1) as a function of the coefficient of restitution a for /i = 1/4, 
uj = 1/2, and (j> = 0.2. The dashed line refers to the first Sonine approximation, the solid line corresponds to the second Sonine 
approximation, while the symbols are the results obtained from Monte Carlo simulations. Here, -D(l) refers to the elastic value 
of the diffusion coefficient consistently obtained in each approximation. 

equations H15() and l|26() and their corresponding Enskog counterparts is the presence of the factors X22 in (|15l) and X12 
in H2t)|l . which are the configurational pair correlation functions for the fluid-fluid and impurity-fluid pairs at contact, 
respectively. A good approximation for these functions is provided by the extended Carnahan-Starling form |46j |: 



X22 



1 - 



(56) 



X12 



1 



2(1 



CT12 2 (1 



- 

\0-i2 



(57) 



where (j) = (7r/6)n2cr| is the solid volume fraction. According to Eqs. (|56|l and H57|l . the correlation functions have 
the same density dependence only when the size ratio is equal to 1. Given that the fluid is in a homogenous state, it 
follows that X22 and X12 are uniform. Thus, it is evident that, when properly scaled, the solutions obtained here for 
the Boltzmann and the Boltzmann-Lorentz equation can be directly translated to the Enskog equation. To illustrate 
the influence of the density on the predictions previously made for a dilute gas, in Fig.|51we plot D/D{1) as a function 
of the coefficient of restitution for = i, a; = i, and a solid volume fraction (p — 0.2. As for a low-density gas, 
the second Sonine approximation clearly improves the results obtained from the first Sonine approximation. We see 
that the second Sonine approximation compares very well with simulation data, except for quite large dissipation 
(a ~ 0.5). 

In summary, there is growing theoretical support for the validity of hydrodynamic transport processes in granular 
fluids. However, some care is needed in translating properties of normal fluids to those with inelastic collisions. As 
seen here, in the case of the diffusion coefficient, the mass transport is strongly influenced by dissipation and this 
dependence may be of some relevance in sedimentation problems, for example. 
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APPENDIX A: EXPLICIT EXPRESSIONS OF d AND a 



In this Appendix, we will quote the expressions of the partial cooling rate Ci and the coefficient ci. These results 
were obtained in Ref. [T]| for arbitrary values of xi. Here, we will display both expressions in the tracer limit xi — > 0. 
The cooling rate Ci is defined by Eq. H22|) . By using the leading Sonine approximations H17|l and H25|l and neglecting 
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nonlinear terms in ci and C2, Ci can be written as [ll| 

Ci = ^10 + Aiici + A12C2, 

where 



AlO = g\/7rn2Cri2WoA*21 



1/2 



M21 



(l + ai2) 1- ^(l + ai2)(l + ( 



(Al) 



(A2) 



All = ^VTrn2al2Vom^^^^ (1 + "12) [2(3 + 46*) - 3/i2i(l + ai2)(l + 9)] , 



(i + ey 



-3/2 



(A3) 



X12 = -Y^VT^n2<7l2VoH2i i—^j (l + ai2)[2 + 3/i2i(l + ai2)(l + 6l)]. (A4) 

Here, 9 — miT /m2Ti — fi/j with /i = mi/m2 and 7 = Ti/T. 
The coefficient ci is defined by 

_ 8 ^ m? /■ 4 ^(0) 15 



dv.Vi<°)-^). (A5) 



15 \4niTi 

This coefficient is detemined by substitution of Eqs. p7|l and H25I) into the Boltzmann-Lorentz equation (|20|1 . multi- 
plying that equation by d*, and integrating over the velocity. When only linear terms in ci and C2 are retained, the 
result is found to be [ll| 

-^^Cifi + ^) =^^lo + ^^llCl + ^!l2C2, (A6) 

where 

ilio = 2^/¥n2gl2^'oA'l2/^2l ' (1 + 012) [-2 (6 + 56*) + /i2i (1 + 012) (1 + 0) (14 + 56 



2 612 ' V 2 



(1 +_^"'/' 

95/2 



-8^x2, (1 + ai2)' (1 + ef + 2fil, (1 + ai2)' (1 + 9f\ , (A7) 

^11 = ^n2(T22T;oAi?2Ai2i-^^^|^(l + «i2) [-2 (90 + 2310 + 18402^40^3) 

+3Ai2i (1 + ai2) (1 + 9) (70 + 1176* + 446|2) - 24/i2i (1 + ai2f (1 + 9f (5 + 46*) 

+3OAiii(l + ai2)^(l + 0)^] , (A8) 

^12 = ^n2(7i^2«oA^i2M2i ^gf/2 ^ + "12) [2 (2 + 50) + 3a^21 (1 + Q12) (1 + 9) {2 + 59) 

-24pili{l + ai2f {l + 9f +3Qfili{l + ai2f {l + 9f] . (A9) 

The final expression of ci is obtained by substitution of Eq. IjAip into Eq. l|A6|l and neglecting nonlinear terms in ci 
and C2. The result is 

_ Aio + A12C2 + ^A^r2^6l2 (»io + ni2C2) 
iAio + All + ^A*r2'^'f^ii 

Once the coefficient ci is given in terms of 7 and the parameters of the mixture, the temperature ratio 7 can be 
explicitly obtained by numerically solving the condition for equal cooling rates: 

Aio + Aiici + A12C2 = C2, (All) 

where C2 is given by Eq. (|19|l . 



16 



APPENDIX B: FIRST AND SECOND SONINE APPROXIMATIONS 

In this Appendix we determine the coefRcients ai and 02 in the first and second Sonine approximation. Substitution 
of Eq. (|39|l into the integral equation (|36|l gives 

- CTdriaifiM^ + a2fi,MSi) - ai J12 [/i,a/v, /2] - 02 Ji2[/i,Af Si, /2] 

Next, we multiply Eq. IBll) by miv and integrate over the velocity. The result is 

{-CTdr + Va)niTiai + niTiV},a2 = -n2T'i- (B2) 
Here, C = C2 is given by Eq. (|19|l and we have introduced the quantities 

mi 



"in 



^j^J d^^r■J,2[hM^.hl (B3) 



mi 
6n 



^ j dvv Ji2[/i,mSi,/2]. (B4) 



From dimensional analysis TiOi ^ T^l"^ so the temperature derivative can be performed in Eq. ljB2p and the result is 

{va - ^C)ai + Vha2 = , (B5) 



where x\ = ni/n2. If only the first Sonine correction is retained (which means 02 ^ 0), the solution to (|B5|I is 



-1 



ai[l] = Ht- (B6) 

J^a- 2C 

Here, ai[l] denotes the first Sonine approximation to ai. Equation (jB6p leads to the expression (|45() for the diffusion 
coefficient D[l] when the second equality in (|42|l is considered. 

To close the problem, one multiplies Eq. 1B2|I by Si(v) and integrates over the velocity. Following identical 
mathematical steps as those made before, one gets 

- CT,-')a, + {ua |C)«2 = -\Afr, (B7) 
2 Zxili 

where ci is given by Eq. HA10|I and we have taken into account that rfa2 ~ T''/^. Moreover, we have introduced the 
quantities 

(B8) 



= y dvSi • Ji2[/i3/Si,/2]. (B9) 



15 niT^ 

In reduced units and by using matrix notation, Eqs. I|B5I) and |IB7|) can be rewritten as 



(BIO) 



2' 

Here, C* = C/t'o, K = '^a/'^o, ^1 = V],/T vq, i/* = Tvc/vq, and v*^ ^ Vd/vo with = ra20-|wo. Further, aj = xiVqQi 
and 02 = xiTi'Qa2. The solution to Eq. ljB10|l provides the explicit expression of the second Sonine approximation 
al[2] to a^: 

n*m-n*U] i<-lC)K-IC-{ci/2jX] 

^ - '^ \K-'2Cm-IC)-K[-*-iC/i)Y ^ ^ 

Equation IjBlip yields directly the expression (|44ll for the second Sonine approximation D[2]. 
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APPENDIX C: EVALUATION OF THE COLLISION INTEGRALS 



In this Appendix we evaluate the quantities Va, i^b, i^c, and Vd defined by the collision integrals l|B3|) . l|B4p . (jBSj), 
and (jB9|) . respectively. Three of them, Va, I'c, and i^d: were already determined in Ref. |25j | for arbitrary composition. 
For the sake of completeness, we display now their explicit expressions in the tracer limit (cci — > 0). In reduced units, 
they are given by 



4 ^ f C7l2 
3 V (72 



(1 + au) 



1 



1/2 



1 



C2 

16 



(CI) 



..=^V^^(^)'(l + ai2/ 

15 ^12 V 



4 ^JJ.li I (Ti2 



1 



1/2 



(C2) 



where 
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1 



-A, 



A, = 5(1 + 2/3) + M2i(l + 0) [5(1 - ai2) - 2(7ai2 - 11)136-^] + IS^^^-^ 
+2^2^(2^22 - 3ai2 + 4)0-1(1 + Of ~ 56'-i(l + 0) 

c ^ 

+ 77 TT-^ {361^21(1 + ai2) [4m2i(1 + ^12) - 5] + 61 [2/ii2 (7m2i(1 + "12) - 5) 
16 (1 + ay 

+M21 (-5(9 + 7ai2) + /i2i(38 + 62ai2 + 24a2^))] - 15 + h^pS.^ ~ 20^21(3 + 0^12) 
+2//2^(40 + 19ai2 + 6a?2) + 2a*12 [M2i(61 + 7ai2 - 20)]} , 



(C3) 



(C4) 
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2^21 



{2oi\2 - 3ai2 + 4) (8 + 56*) 



-A*2i(l + Q) [2/361-2(8 + 56')(7ai2 - 11) + 26i-i(29ai2 - 37) - 25(1 - ^12)] 
+ 18/3261-2(8 + 561) + 2/361-1(25 + 666*) + 56'-i(6 + II6I) - 5(1 + 6')0-2(6 + 561) 

+ ^ (1 + (?)-' {150^21(1 + ai2)(4M2i(l + ai2) - 5) 
16 

+2 [45 + 540AiL + 16M2i(ai2 - 36) + 4^2^(134 + 5ai2 + 60^2) 
-4^12(148 + Ai2i(7ai2 - 263))] + 0^ [-30 - A*2i(267 + 217ai2) 
+ 14/^2^(17 + 29ai2 + 12a?2) + 10/^12 (7/i2i (1 + ^12) - 5)] 

+61 [-315 + 270/1^2 - 2M2i(55ai2 + 57) + Ai2i(440 + 326ai2 + 156a^2) 
+2/ii2(-2 + ^2i(7ai2 + 277))]}. 



(C5) 



In the above expressions, /5 — [i\2 ~ A'216' and C2 and are given by Eqs. ((TH|l and respectively. 

It only remains to evaluate I't, which is defined by the collision integral (|B4p . To simplify the integral, a useful 
identity for an arbitrary function h(y\) is given by 



y"dvi/l(vi)Ji2 [vil/1,/2] = (Ti^ / ^'^l y C?V2/l(vi)/2(v2) 

X y d?e(^-g)(a-g) [/i(v'/)-/i(vi)], 



with 



v" = Vi - /i2l(l + Q!12)(5' • g)CT , 

and g = vi — V2 is the relative velocity. Use of (|C6p in Eq. IIB4|) gives 



^^^^^12^21(1 + 0^12) y" rfvi y dv2g/i,Af(vi)/2(v2) [Si(vi) • g] . 



(C6) 



(C7) 



(C8) 
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Substitution of the distribution function /2 from Eqs. H16|) and H17|l gives 



37r2 



Ml2 



^12 
0-2 



C2 



*4 r *2 
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(C9) 



where v* = v/vq and g* — g/wo- The integrals appearing in (|C9|) can be evaluated by the change of variables 

x = v*-v;, y = 0vt+v;, (CIO) 



with the Jacobian (1 + ^ 
quantity = Vb/Tvo is 



The integrals can be easily performed and the final expression for the dimensionless 
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